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A  correlation  analysis  of  a  thickness  tendency  equation  including 
heating  and  vertical  motion  terms  was  carried  out.   One  of  the  main  pur- 
poses of  the  investigation  is  to  provide  a  firmer  foundation  for  a  1000- 
500  mb  thickness  forecast  scheme. 

A  large  data  sample  was  used  and  single  and  multiple  linear  coeffi- 
cients were  calculated.   These  were  computed  for  a  wintertime  situation 
with  separate  analyses  for  Atlantic  and  Pacific  areas. 

The  resulting  correlation  coefficients  for  the  advection,  vertical 
motion  and  heating  terms  gave  substance  to  the  hypothesis  that  the  observed 
change  in  thickness  may  be  fairly  well  approximated  by  the  change  due  to 
horizontal  advection  only.   The  vertical  motion  term  gave  a  smaller  but 
significant  correlation,  while  heat  exchange  terms  did  not  correlate 
significantly. 

The  author  wishes  to  express  his  appreciation  to  Professor  G.  J. 
Haltiner  for  his  guidance  and  encouragement,  to  Mrs.  Jean  Bow  of  the 
computer  center  at  the  U.  S.  Naval  Postgraduate  School,  for  her  help  in 
Fortran  programming,  and  especially  to  Mr.  James  R.  Clark  of  Meteorology 
International,  Inc.,  for  his  assistance  in  data  collection  and  map  re- 
production.  The  author  also  wishes  to  express  appreciation  to  the  U.  S. 
Naval  Fleet  Numerical  Weather  Facility  for  making  their  data  and  computer 
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e  Water  vapor  pressure  of  air  (mb) 
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e  Saturated  water  vapor  pressure  at  the  temperature 

of  the  water  surface  (mb) 

h  1000-500  mb  thickness 

h  Combined  local  change  of  thickness  and  estimated 

horizontal  advection 
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1000-500  mb  thickness  tendency 
H  Total  heat  exchange  over  12-hour  period 
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1 .    INTRODUCTION 

The  ability  to  produce  analyses  and  prognoses  of  weather  maps  by 
numerical  methods  using  high  speed  digital  computers  on  an  operational 
basis  is  now  a  reality.  Most  of  the  efforts  in  numerical  weather  pre- 
diction have  been  in  relation  to  the  pressure  and  temperature  fields; 
however  some  pioneering  attempts  have  been  made  to  predict  precipitation, 
clouds  and  visibility. 

Perhaps  the  most  successful  forecasts  have  been  of  the  500-mb 
contours.   However,  somewhat  less  success  has  been  achieved  in  obtaining 
a  good  surface  or  1000-mb  prognosis.  Many  schemes  have  been  tried  with 
varying  degrees  of  success.   One  method  that  has  a  great  deal  of  promise 
makes  use  of  the  500-mb  and  1000-500  mb  thickness  prognoses  to  obtain  a 
1000-mb  forecast.   The  U.  S.  Naval  Fleet  Numerical  Weather  Facility 
(FNWF)  has  utilized  an  empirically  adjusted  advective  technique  for 
thickness  prediction.   A  more  complete  thickness  forecast  model  includ- 
ing vertical  motion,  heating  and  friction  is  presently  under  investiga- 
tion by  Professor  G.  J.  Haltiner  of  the  U.  S.  Naval  Postgraduate  School. 
This  program  utilizes  a  thickness  tendency  derived  from  the  first  law  of 
thermodynamics  and  is  shown  in  a  subsequent  section  of  this  paper.   In 
some  of  the  earlier  applications  of  this  equation  some  of  the  terms  have 
been  omitted  by  way  of  simplification  because  of  the  lack  of  information 
necessary  to  calculate  these  terms,  for  example,  diabatic  heating.   Re- 
cently FNWF  has  made  computations  of  the  heat  exchange   on  a  regular 
basis.   With  this  data  now  available  it  is  appropriate  to  consider 
whether  a  better  estimate  of  a  thickness  tendency  might  be  obtained  if 
more  terms  of  the  equation  were  retained  and  fewer  simplifying  assumptions 


made. 

The  purpose  of  this  investigation  is  to  conduct  a  first  step  in  a 
determination  of  the  usefulness  of  this  equation,  that  is,  to  statisti- 
cally correlate  the  observed  1000-500  mb  thickness  changes  with  estimates 
of  the  various  terms  in  the  thermodynamic  equation.   As  shown  by  Panofsky 
and  Brier  [8]  linear  correlation  coefficients  give  a  measure  of  associa- 
tion between  two  variables  which  does  not  depend  upon  any  arbitrary  choice 
of  units  by  which  the  original  variables  were  measured.   From  the  correla- 
tion coefficients  obtained  here  each  term  will  be  examined  as  to  its  re- 
lation to  the  observed  thickness  change  and  its  significance  as  a  pre- 
dictor.  From  these  results  the  value  of  this  equation  in  a  prediction 
scheme  can  be  estimated. 
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2.   AREA  OF  INVESTIGATION 

FNWF  routinely  analyzes  data  over  the  Northern  Hemisphere  on  a 
63  x  63  grid,  with  a  grid  distance  of  381  KM. 

As  shown  later  in  this  paper  the  scheme  used  in  obtaining  upwind 
points  for  data  extraction  was  designed  for  use  on  a  22  x  22  grid  extract 
ed  from  any  section  of  FNWF's  63  x  63  grid.   The  total  heat  exchange 
calculations  were  originally  derived  for  use  over  the  oceans  only,  how- 
ever FNWF  calculates  this  term  over  the  whole  grid,  and  then  plots  only 
the  ocean  points.   Because  the  land  area  points  were  to  be  kept  at  a 
minimum,  all  available  values  of  the  heat  exchange  were  used.   The  known 
decrease  in  the  accuracy  of  the  value  over  land  as  compared  to  that  over 
the  ocean  nevertheless  was  accepted  in  this  study  and  recognized  as  a 
possible  source  of  error. 

Due  to  the  facts  mentioned  above  two  different  areas  of  22  x  22  grid 
size  were  selected  for  this  study.   These  areas  were  chosen  because  most 
of  the  grid  is  oceanic  with  a  minimum  of  continental  area. 

Area  1,  which  will  be  referred  to  as  the  Atlantic  area,  is  shown  in 
figure  1.   This  area  includes  the  North  Pole  and  the  Gulf  Stream  and  has 
a  minimum  of  equatorial  area;  the  latter  is  near  the  boundary  of  FNWF  's 
grid  where  results  are  considered  to  be  less  accurate. 

Area  2,  which  will  be  referred  to  as  the  Pacific  area,  is  shown  in 
figure  2.   This  area  is  almost  entirely  oceanic,  but  does  contain  more 
of  the  equatorial  or  boundary  area.   These  areas  will  be  treated  separ- 
ately and  the  correlation  coefficients  will  be  calculated  independently. 
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3 .    BACKGROUND 

The  thickness  prediction  equation  may  be  derived  in  the  following 
way.  The  first  law  of  thermodynamics  can  be  written 

i        ii        iii  iv       v  (1) 

Combining  the  equation  of  state  and  the  hydrostatic  equation  gives 


T*-  tf 


37  TJ*  (2) 


Substituting  (2)  into  (1)  in  terms  (ii)  and  (iii)  gives 
■? 


±,.^in)-^ v.v(if)  +  ^(^-§)    o) 


Rearranging,  we  have 


ttr  -   v  v^rp  +  <*t  W  Cft)    c^+  (4) 


Integrating  between  1000  and  500  mb 


j^-JV-V^^f/^aJ^     ~ff&^  (5) 


^  '/<POO 


where 


r-.- A-(ft-S5.) 


Now  integrating  over  a  12-hour  time  period  gives 


Jhsoo        .  fiasco  /'V*6« 


12 


The  purpose  of  this  investigation  is  to  correlate  the  12-hour  thick- 
ness changes  with  numerical  approximations  of  the  terms  constituting  the 
right  hand  side  of  equation  7.   Thus  Ah.      will  be  correlated  with 
estimated  mean  values  of  the  horizontal  advection,  A,  the  vertical  motion 
term,  VJ,  and  the  total  heat  exchange,  H,  for  the  12-hour  period.   This 
relationship  may  be  written 


Ak   =    A    +    W    +   H 


(8) 


The  coefficients  of  the  terms  in  equation  7  can  be  considered  as 
merely  constants  of  proportionality  and  thus  can  be  omitted  in  the 
statistical  evaluation. 

The  parameter   0~   as  defined  here  differs  slightly  from  the  value 
of   0""   used  by  FNWF.   However,  here  again  the  difference  can  be  con- 
sidered essentially  as  a  constant  of  proportionality  and  will  not  affect 
this  study. 
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4.    DATA  COLLECTION 

The  basic  data  for  this  study  was  obtained  for  the  three-day  period 
9  to  11  February  1966  from  the  Fleet  Numerical  Weather  Facility. 

The  initial  fields  utilized  were  the  SD-500  field  which  is  the  dis- 
turbance component  (perturbation)  of  the  500-mb  height  field,  the  SR-500 
field  which  is  the  residual  component  (basic  flow)  of  the  500-mb  height 
field,  the  SD-1000  field  which  is  the  disturbance  component  at  the  1000 
mb  level  and  the  SR-1000  field  which  is  the  basic  flow  at  the  1000-mb 
level. 

These  fields  were  used  in  a  routine  originally  designed  for  the 
interpretation  of  meteorological  satellite  observations  as  developed  by 
Nagle,  Clark  and  Holl  [7]  and  modified  by  Clark  for  use  in  obtaining  the 
advective  change  in  thickness,  A.   The  scheme  entails  generation  of 
stream  function  fields,  ^  ,  figures  3  and  4,  from  the  SD  and  SR  fields 
averaged  over  the  1000-500  mb  layer.   By  a  timewise  interpolation  of  the 
12-hour  distributions,  hourly  distributions  of  r  were  obtained  which 
would  act  as  evolving  frame-  of-  reference  steering  fields.   From  these 
T  fields  a  horizontal  velocity  distribution  was  computed.   This  wind 
was  based  on  the  geostophic  approximation  and  in  this  case  was  the  mean 
wind  in  the  layer. 

At  this  time  the  12-hourly  winds  for  a  22  x  22  subsection  of  the 
63  x  63  grid  were  extracted.   This  sub region  was  then  interpolated  in 
both  time  and  space  to  give  hourly  winds  into  an  expanded  44  x  44  array 
based  on  the  22  x  22  section.   This  maneuver  quadruples  the  number  of 
data  points  contained  in  the  region. 

By  utilizing  the  stream  function  and  the  wind  field,  an  upstream 
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trajectory  was  projected  and  the  position  of  a  column  of  air  at  time 

t  -  12  was  determined.  This  was  done  for  each  data  point  in  the  region. 

Checks  in  the  program  flagged  any  point  whose  upstream  position  was  out 

of  the  44  x  44  grid.  This  point  was  then  discarded  for  correlation 

purposes. 

Several  steps  were  involved  in  the  procedure  used  to  determine  A. 
Thickness  fields,  figures  5  and  6,  were  obtained  from  FNWF's  D-500  and 
D  1000  fields.   These  D  fields  give  the  difference  between  the  height 
of  the  surface  (500  and  1000  mb,  respectively)  and  its  corresponding 
standard  height.  The  trajectory  scheme  was  then  applied  to  this  field 
to  yield  a  new  parameter,  h,  which  is  the  combined  local  change  of  thick- 
ness, Alr\,   ,  and  the  estimated  horizontal  advection,  A.   Figures  7  and 
8  depict  this  parameter  in  the  Atlantic  and  Pacific  areas,  respectively. 
The  strong  gradient  in  the  extreme  western  section  seen  in  these  and 
subsequent  maps  was  due  to  the  check-flag  system  discussed  previously. 
These  are  the  areas  of  inflow  where  the  upstream  position  was  off  the  grid 
causing  a  zero  to  be  inserted.  While  the  zero  value  was  eliminated  be- 
fore correlations  were  made,  the  point  was  plotted  as  such  resulting  in 
tight  gradients  where  none  actually  exist. 

The  local  or  observed  12-hour  change  in  thickness,  figures  9  and 
10,  was  produced  by  subtracting  the  field  at  time  t  -  12  from  the  field 
at  time  t  .  To  get  an  estimate  of  A  at  each  point,  Ah,    was  subtracted 
from  h;  this  was  not  completed  until  the  format  of  the  data  had  been 
changed  and  the  correlation  program  was  in  use. 

A  further  modification  of  the  routine  allowed  W  to  be  computed. 
This  vertical  motion  term  was  obtained  by  taking  the  product  of  0~  and 
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Uj   ,  the  vertical  motion  in  (x,y,p,t)  coordinates,  and  averaging  in 
both  time  and  space. 

To  accomplish  this  sigma-1  and  sigma-3,  the  static  stability  for  the 
layers,  1000-750  and  600-750  mb  respectively,  weie  utilized.  These  0~  's 
were  derived  from  a  model  by  Hall,  Bibbo  and  Clark  [2]  which  is  charact- 
erized by  pressure  levels  defining  a  number  of  layers  within  which  the 
static  stability  is  independent  of  pressure. 

These  terms  were  multiplied  by  the  U)     fields  at  850  and  500  mb. 
FNWF  computes  these  parameters  by  a  numerical  solution  of  a  quasi- 
geos trophic  omega  equation,  as  shown  by  Haltiner,  Clarke  and  Lawniczak 
[1].   The  resulting  fields  were  0~oJOCft,  figures  11  and  12,  and  <TiO        , 
figures  13  and  14.   The  trajectory  scheme  was  again  used  to  obtain  up- 
wind values  of  these  terms  at  time  t  -  12  which  were  then  averaged  with 

o 

values  at  t  to  yield  the  12-hour  averages,  <TUJ  and  <3"CJcrir>,  shown 

O  ojU         5UU 

in  figures  15,  16,  17  and  18.   The  final  averaging  of  <rtUft_n  and  G~U)   -- 
to  obtain  W  was  completed  by  the  correlation  program  itself. 

The  same  procedure  was  used  to  obtain  H.   This  term  was  derived  from 
the  fields  of  Q  ,  the  total  air-ocean  interchange  as  produced  by  FNWF 
shown  in  figures  19  and  20.   The  computation  program  of  the  total  heat 
exchange  by  numerical  methods  and  its  adaption  to  the  computer  evolved 
from  earlier  manual  methods  developed  by  Laevastu  [4].  The  heating  func- 
tion was  obtained  from  the  following  equation: 

Qru  -    QS  -  Gr    -  Qb   -   Qk   -  Qe  (9) 

-2        -1 
where  all  Q's  will  be  expressed  in  units  of  g-cal  cm   (24  hr.) 

Because  of  the  relative  newness  of  the  computerized  form  of  this  product 
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each  term  will  be  briefly  discussed. 

The  insolation,  Qs,  is  calculated  from  the  equation 

Qs  =  o.ih-  AK-td  (  \  -  o.ooo^C  )  (io) 

which  takes  into  account  the  total  time  that  the  sun  is  above  the  horizon, 

t,,  the  noon  altitude  of  the  sun,  A  ,  and  a  cloudiness  factor,  C.  Since 
d  n 

it  is  virtually  impossible  to  take  into  consideration  various  types  of 
clouds  and  their  different  effects  on  the  insolation,  Laevastu  [3] 
devised  an  average  cloudiness  factor  which  requires  only  a  knowledge  of 
total  cloud  tenths.   The  results  give  a  reasonable  estimate  of  the  total 
incoming  radiation. 

The  albedo,  Q  ,  is  an  estimate  of  the  radiation  reflected  from  the 
surface.  Laevastu  [5]  statistically  determined  an  empirical  formula  for 
calculating  this  term. 

Qr  =  0.15  Qs  -  (o.o\  Qs)2  (11) 

This  formula  is  valid  only  for  daily  calculations.   A  simplified  formula 
for  the  computation  of  short  term  (hourly  or  3-hourly)  albedo  is 

Qr  =  Q5  -=r-  ,    °<-    =  Average  solar  altitude         (12) 
The  effective  back  radiation,  Q,  ,  is  the  long  wave  radiation  given 
off  by  the  sea  surface.  This  represents  a  heat  loss  at  the  surface 

Qb  *  (^7  -  l.*fe  7^  -  0.9i-O)(  I  -O.OlfcS-C)     (13) 

and  makes  use  of  a  linear  formula  of  Lonnquist  f6]  which  is  corrected  by 
Moeller's  cloudiness  factor. 
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The  convective  transfer  of  sensible  heat,  Q,  ,  was  obtained  from 

h 

Q^  =  3^(o.2fc+o.m  v)(tw  -To.)   ;  7Z   >  Tl  (14) 

when  heat  is  transferred  from  the  water  to  the  air;  and 

Qv^r  3V(T„  -Tou)  }    Tw  <  T«_  (15) 

when  heat  is  transferred  from  the  air  to  the  water. 

The  latent  heat  transfer,  Q  ,  is  that  heat  addition  or  loss  due  to 

e 

evaporation  or  condensation.  Lavastu  [5]  shows  that  a  modified  formula 
of  Rohwer  has  been  found  to  be  most  accurate  in  estimating  evaporation; 
namely, 

A  modified  form  of  this  equation  can  be  used  for  condensation  at  the  sea 
surface, 

Qc=    0.077  v  (o,<}ee„  -  eA)  Lt  p  rf  o,i*e»<^  (i7) 

To  calculate  these  terms  and  get  a  total  heat  exchange  product 
certain  basic  synoptic  numerical  analyses  must  be  available.   These  are: 
Surface  winds,  surface  air  temperature  and  water  vapor  pressure  and 
clouds.   These  are  all  routinely  produced  at  FNWF. 

By  utilizing  the  trajectory  program  again,  the  t   -  12  values  of 
Q  were  obtained.  Averaging  over  the  12  hours  produced  H,  figures  21 
and  22.   No  further  data  was  needed. 
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5.    OBJECT  AND  METHOD  INVESTIGATION 

The  purpose  of  this  investigation  was  two-fold.   An  attempt  was 
made  a)  to  statistically  correlate  the  observed  change  in  thickness 
with  estimated  values  of  the  terms  on  the  right  hand  side  of  equation 
(7)  and  b)  to  estimate  the  significance  of  each  term  in  a  prediction 
scheme. 

To  accomplish  this  a  multiple  regression  and  correlation  analysis, 
the  BMD03R  programmed  for  the  CDC  1604  electronic  computer  and  converted 
for  COOP  MONITER  Fortran  63,  was  utilized.   This  program  will  perform 
linear-multiple  correlation  analysis  on  data  in  any  combination  by  allow- 
ing the  designation  of  a  dependent  variable  and  deletion  of  sets  of  in- 
dependent variables.   This  can  be  done  as  many  times  as  desired. 

At  the  completion  of  the  data  collection  phase  there  were  five 
parameters  stored  for  each  data  point  over  a  total  of  8229  Atlantic  and 
6739  Pacific  points  which  were  considered  good  for  correlation  purposes. 
These  variables  were:  (1)  £  ;  (2)^Jg50;  O)™*  500'>    (4)  H  ;  and  (5)  .Ak  . 
To  place  this  data  in  a  form  appropriate  for  the  BMD  program  the  format 
had  to  be  changed  from  a  row-wise  listing  to  a  column  presentation. 
After  accomplishing  this,  spot  checks  were  made  to  ensure  that  the  data 
had  not  gotten  out  of  phase  during  the  transformation. 

Since  the  parameters  A  and  W  had  not  actually  been  obtained  as  yet, 
the  transgeneration  feature  of  BMD03R  was  used.   This  allows  for  the 
generation  of  new  variables  by  performing  certain  mathematical  opera- 
tions on  the  original  variables.   Here  A,  variable  6,  was  obtained  by 
the  subtraction  (1)  -  (5)  =  (6)  and  W,  variable  7,  was  produced  from 
[(2)  +  (3)]%  =  (7).   The  data  was  then  ready  for  the  correlation  analy- 
sis. 
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6.    RESULTS  AND  DISCUSSION 

The  data  was  divided  into  Atlantic  and  Pacific  regions  with  identi- 
cal calculations  made  for  each  area.   Initially,  correlation  analyses 
were  carried  out  to  show  relationships  between  pairs  of  variates.  Many 
combinations  of  independent  and  dependent  variables  were  used.   The  re- 
sults are  shown  in  Tables  1  and  2. 

After  examining  the  figures,  multiple  linear  correlations  were  made. 
This  was  done  by  progressively  adding  the  term  from  equation  (8)  with  the 
highest  partial  correlation  coefficient. 

As  a  side  operation,  a  multiple  correlation  analysis  was  made  on 
the  terms  of  the  equation 

h\    <     W  +   ft 

These  results  are  shown  in  Tables  3  and  4. 

The  values  of  R  obtained  from  the  Atlantic  data  were  slightly  higher 
than  those  in  the  Pacific  in  almost  every  case.   This  can  be  partially 
explained  by  a)  the  scarcity  of  reporting  stations  in  the  Pacific  result- 
ing in  larger  interpolated  areas,  and  b)  the  Pacific  area  contains  a 
larger  band  of  equatorial  or  boundary  region  where  the  analysis  is  known 
to  be  less  accurate. 

The  coefficient  obtained  in  the  linear  correlation  of  ^K.  and  H  was 
extremely  low  in  both  the  Atlantic  and  Pacific. 

The  coefficient  resulting  from  the  correlation  of  Ah, and   W  was  lower 
than  expected  but  was  approximately  the  same  in  both  Atlantic  and  Pacific. 
^  h*   was  correlated  with  CW    and  <fU)  _nn  separately  to  see  if  it  was 
necessary  to  obtain  W  or  if  the  value  of  the  parameter  obtained  at  a 
single  level  would  be  sufficient.   In  the  Atlantic  there  is  little 

20 


2 
Dependent  Independent  R  R 

Variable  Variable 

Ak  "S  -8599         .7394 

AK.  W  .1936         .0375 

A^  H  .0655         .0043 


AK^  ^^aso  -1729       -0299 


AK,  C=&Jroo  -1850         .0342 

W  A  -4513         .2037 

"W  TT  .2200         .0484 


A  H  -1799         .0324 

\\,  ~\4  .2425         .0588 

Irv.  W  -5802         .3366 


Table  1.  Results  of  non-multiple  linear  correlation  analyses, 
Atlantic.   R  is  the  correlation  coefficient  and  R2  is  the  co- 
efficient of  determination. 
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Dependent 
Variable 

Ak 
Ah- 

AW, 
AH- 


Independent 
Variable 


w 

Vl 


A 
W 


croj 


?50 


"A 


H 
W 


7967 


1752 


0985 


2440 


0825 


2569 


0404 


0962 


0483 


5373 


.6347 
.0302 
.0097 
.0595 
.0068 
.0660 
.0016 
.0093 
.0023 
.2886 


Table  2.   Results  of  non-multiple  linear  correlation  analyses, 
Pacific. 
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2 
Dependent  Independent  R  R 

Variable  Variable 

AW,  "A  .8599         .7394 


AW,  A  +  VA/  .8871         .7869 


AW.             A,  +  W  +  W                           8889  #7901 

ft.                  "W+~\T              .5920  .3505 

Table  3.  Results  of  the  multiple-linear  correlation  analysis, 
Atlantic. 

2 

Dependent              Independent             R  R 

Variable  Variable 

AW.                   A                .7967  .6347 

AH-                A  +W             .8562  .7331 

AK              X^W-i-"ff           .8625  .7439 

|^                  W-rlT             .5418  .2936 

Table  4.  Results  of  the  multiple-linear  correlation  analysis, 
Pacific 
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difference  in  the  coefficients  but  in  the  Pacific,  the  three  differ 
considerably. 

The  multiple  correlation  coefficient  in  the  Atlantic  increased 
from  .8599  to  .8871  to  .8889  as  W  and  H  were  added  to  "A.   The  increase 
in  the  Pacific  was  somewhat  greater,  going  from  .7967  to  .8562  to  .8625 
as  these  terms  were  added.   Here  the  inclusion  of  the  vertical  motion  is 
more  significant. 

The  analysis  of  h>  with  W  and  H  was  attempted  with  the  thought 
that  K.  might  be  used  in  a  cloud  forecast  scheme  currently  under  study. 
Although  the  correlations  in  both  areas  are  significant,  they  are  certain- 
ly not  high  enough  to  consider  replacing  vertical  velocities  by  H-  in  a 
cloud  forecasting  scheme. 

The  remaining  analyses  clearly  show  no  significance  and  will  not  be 
discussed  further. 
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7.    CONCLUSIONS  AND  RECOMMENDATIONS 

As  a  result  of  this  study,  which  was  made  using  wintertime  data 
only,  the  following  conclusions  were  reached. 

(1)  Observed  changes  in  thickness  over  a  12-hour  period  are  due 
almost  entirely  to  horizontal  advection,  at  least  as  approximated  in  this 
study.   This  result  lends  strong  support  to  the  simple  advective  scheme 
for  predicting  thickness  which  has  been  used  operationally  by  FNWF  dur- 
ing the  past  several  years,  i.e. 

Ah.    *c  V-  7  h, 

(2)  The  vertical  motion  term,  W  ,  though  of  secondary  importance, 
does  contribute  significantly  to  the  change  in  thickness  and  should  be 
retained  in  a  prediction  scheme  for  increased  accuracy.   Of  course,  a 
decision  must  be  made  as  to  whether  the  improvement  in  forecasting  skill 
is  worth  the  time  and  effort  necessary  to  include  vertical  motion. 

(3)  The  heat  exchange  term,  H,  as  presently  being  calculated  at 
FNWF  and  as  utilized  in  this  study  contributes  virtually  nothing  to  the 
change  in  thickness  and  should  be  omitted  in  a  thickness  prediction  scheme, 

As  pointed  out  earlier,  this  study  dealt  only  with  winter  d^i:,;.   An 
obvious  and  necessary  extension  should  be  the  application  to  data  from 
the  other  seasons,  especially  summer. 

It  is  hoped  that  the  results  of  this  study  will  provide  guidance  to 
individuals  or  groups  engaged  in  the  development  of  thickness  prediction 
schemes. 
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